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ABSTRACT 

We have tested several models of GRB luminosity and redshift distribution 
functions for compatibility with the BATSE 4B number versus peak flux 
relation. Our results disagree with recent claims that current GRB observations 
can be used to strongly constrain the cosmic star formation history. Instead, 
we find that relaxing the assumption that GRBs are standard candles renders 
a very broad range of models consistent with the BATSE number-flux relation. 
We explicitly construct two sample distributions, one tracing the star formation 
history and one with a constant comoving density. We show that both 
distributions are compatible with the observed fluxes and redshifts of the bursts 
GRB970508, GRB971214, and GRB980703, and we discuss the measurements 
required to distinguish the two models. 

Subject headings: cosmology: observations - gamma rays:bursts - stars:neutron 
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1. Introduction 

The recent measurements of redshifted absorption lines in the optical counterpart 
associated with GRB970508 ( [Metzger et al. 1997|) , and emission lines from the galaxies 



associated with GRB971214 QKulkarni et al. 199SQ and GRB980703 ( P.jorgovski et al. 1998Q 



have confirmed the cosmological origin of gamma ray bursts. Such measurements offer a 
new method of estimating the distance and luminosity distribution of GRBs. Determining 
the redshift distribution, in particular, is of enormous importance in distinguishing between 
GRB source models. Much attention has focused on models in which burst activity is tied 
to the massive star formation rate, either directly (e.g., the "failed supernova" models of 
Woosley 1993), or with a time delay (e.g., the in-spiraling neutron star model of Pacyhski 
1986). In such models, GRB activity should trace the cosmic star formation history ( Madau| 
19971 ; |b'all, Gharlot, fc Pei 199(j| ). 



Previous analyses have attempted to constrain the GRB distribution using the 
cumulative number versus peak flux, iV(< S), relation observed by the BATSE instrument 
(e.g. [Rutledge, Hui, fc Lewin 1995| ). By making additional assumptions about the GRB 
luminosity distribution (usually that bursts are "standard candles"), various authors have 
reached strong conclusions about the GRB distribution. For example, Totani (1998) recently 
used GRB observations to argue that our understanding of the cosmic star formation rate 
is incorrect. 

As we show in §|3|, relaxing the constraints on the intrinsic GRB luminosity function 
greatly reduces the power of GRB number counts to constrain the burst distribution. 
Standard candle models are no longer tenable; there are now redshift measurements of three 
bursts that vary in intrinsic luminosity by a factor of 20.Q We demonstrate that a wide 
range of burst distributions are allowed by the current data, including both the BATSE 
measurements and the measured burst redshifts. However, we show in §^ that further 
redshift measurements, as well as deeper number counts, can be used to "deconvolve" the 
unknown redshift and luminosity distributions. 



2. Theory 

2.1. Tests Using Number Versus Flux 

A source at redshift z has an integrated comoving luminosity in the bandpass 
at the source given by L = ^rdv, where dL/du is the spectrum of the source. If we 
assume a spectral form dL/dv = Cv~ a , we can compute the luminosity L' at Earth in the 



Unless otherwise stated, all models assume an H = 75 km s" 1 Mpc" 1 , fi = 1.0, A = 
cosmology, with no beaming. 
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bandpass [u[, v' 2 ], where v' — u/(l + z), using 

L = / Cv- a dv' = (1 + z)- a / Cz/- Q c/z/ = (1 + z)- a ■ L. (1) 

1 + Z Ju[ Ju! 

The additional factor of 1 + z in equation ([!]) comes from time dilation at the burst source. 
The flux from the source observed at Earth is ([Rutledge, Hui, fc Lewin 1995|) 

S(zL) L(l + z)- a 

The coordinate distance, r(z), is given by Weinberg (1972) 

zqo + (<?o - 1) (y/2q z + 1 
H R q 2 {l + z 

Here i?o, -f^o, and go are the present-epoch expansion, Hubble, and deceleration parameters. 



r ( z ) = u D In I ^ -■ (3) 



A model for GRB sources consists of a distribution in comoving redshift space, dp/dz, 
and a distribution in luminosity space, dg/dL. The redshift distribution has three parts: a 
geometric factor that accounts for the differential volume available at redshift z, a factor to 
include relativistic effects, and a factor n(z) for the comoving differential density of GRBs 
per unit volume at redshift z. Thus, dp/dz is given by 

dp dV 1 . , r(z) 2 dr(z)n(z) , . 

— — ■ = — n(z)=Au; . = =— (4) 

dz dz 1 + z i _ kr(z) 2 "-Z 1 + z 

Here, dV/dz is the comoving volume element per unit redshift, the factor 1/(1 + 2) comes 
from time dilation at the burst source, and k is —1, 0, or 1 for open, flat, and closed 
cosmologies, respectively. 

To test a model, we generate a large number of bursts. We select redshifts from the 
distribution dp/dz, and luminosities from dg/dL. Using equation (0), we determine the flux 
for each burst in the sample, correct for flux- dependent instrumental sensitivity, and form 
a cumulative distribution function iV(< S). We compare the distribution function for the 
model to the observed iV(< S) relation using a Kolmogorov-Smirnov (KS) test, calculating 
a confidence level at which that model can be ruled out. 



2.2. Tests Using Redshift Measurements 

With redshift measurements, we can compare a model's predicted redshift distribution 
to the observed distribution. It is necessary to correct the predicted distribution for the bias 
introduced by a finite flux cutoff. To do this, we generate a large number of bursts from 
the model, correct the sample for instrumental sensitivity, and simultaneously compare 
the redshift and flux distributions using a two-sided, two-dimensional KS test ( press et~al. 
19921) . 
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3. Maps of Parameter Space 

We examine range-dominated power law luminosity distributions, which have previously 
been studied by Emslie & Horack (1994) , Hakkila et al. (1996) , Horack et al. (1996) , 
Horvath et al. (1996) , and Loredo & Wasserman (1998) . These distributions have the form 
dg/dL = L p in the range [L m i n , L max ] and are zero everywhere else. For the star formation 
rate, we use the closed-box model of 280 nm cosmic emissivity, including dust obscuration, 



taken from Fall et al. (1996). This paper and others (e.g. |Madau 1997|) include several 
other estimates of the cosmic star formation history. We choose this model as an example; 
other SFR models would give numerically different but qualitatively similar results. 

First, we find the comoving luminosity of the three GRBs in the BATSE bandpass 
(50-300 keV). GRB971214 has a peak flux of 1.9538 photons cm~ 2 s" 1 on the 1024 ms 



time-scale, and a measured redshift of 3.42 QKulkarni et al. 1998| ). GRB970508 has a peak 



flux of 0.969 photons cm 2 s 1 , absorption features at z = 0.835 ( [Metzger et al. 1997] ) , and 



it most likely in a galaxy at that distance QPian et al. 1998] ). GRB980703 has a peak flux 



of 2.42 photons cm 2 s 1 , and is associated with an object at z — 0.967( [D]orgovski et al, 
I998|) . We obtained the BATSE peak fluxes from Meegan et al. (1998 



To use equation (pD, we need to know a, the GRB spectral index. GRB spectra vary 
considerably ( [Band et al. 1993]) ; a thorough treatment would assume a distribution of 
underlying spectra in its analysis. However, as we will see, the current data do not yet 
justify this expansion of the model phase space. Instead, we simplify by using as our 
single spectral index the recent estimate a = 1.1 ([Mallozzi, Pendleton, fc Paciesas 1996j ). 
With this assumption, we calculate luminosity Li = 2.15 x 10 50 erg/s for GRB970508, 
L 2 = 4.60 x 10 51 erg/s for GRB971214, and L 3 = 6.70 x 10 50 erg/s for GRB980703. 

We examine the section of parameter space Li/10 < L min < Li, L 2 < L max < L 2 ■ 10, 
— 1 < p < —2.5, on a 21 by 21 by 11 grid evenly spaced in the logarithm of L min and L max , 
and evenly spaced in p. 

The BATSE GRB catalog ([Meegan et al. 1998|) , as of late July, includes 2195 bursts. 
To insure uniform efficiency, we restrict our study to the 1277 bursts triggered on channels 
2+3 with a threshold of 5.5 sigma on all three BATSE time scales. We eliminate from this 
list 307 GRBs with peak flux below 0.42 photons cm 2 s -1 , where BATSE's efficieny is near 
unity. Efficiency drops radidly below this cutoff, and uncertainties in efficiency contribute 
significantly to estimated number counts QPendlcton, Hakkila, fc Meegan 1998| ). Using a 
KS test, we compare the remaining 970 bursts to 5000 GRBs generated from each model. 

The resulting confidence levels are shown in Figure [l|. The data clearly show there are 
a number of luminosity distributions consistent with a redshift distribution that goes as the 
star formation history. Our results are consistent with previous work indicating that the 
BATSE data constrain most GRBs to a single decade in luminosity QHorack et al. 19~96| ). 
Despite this constraint, the BATSE data presently allow a very wide range of possibilities. 



- 6- 




Lmin factor Lmin factor 

Fig. 1. — The light contour shows models with significance greater than 0.05, and the dark 
contour shows models with significance greater than 0.33. One minus the significance value 
is the confidence level at which the model can be ruled out. L min factor and L max factor 
indicate the amount by which the luminosities of GRB970508 and GRB971214 are multiplied, 
respectively. Values of p, the exponent of the luminosity distribution, are indicated in each 
panel. 
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4. Case Study of Two Models 

To determine what measurements can differentiate the models, we study two models, 
which we call the star formation rate model (SFR) and the constant comoving density 
model (CCD). SFR uses the star formation rate as its redshift distribution, and has 
Lrmn = 1-52 x 10 50 erg/s, L max = 1.16 x 10 52 erg/s, p = —1.9 for its luminosity distribution. 
The CCD model has a constant comoving redshift distribution, an L m i n = 1.13 x 10 50 erg/s, 
I max = 1.14 x 10 52 erg/s, p = —2.30 for its luminosity distribution. Both these models fit 
the BATSE data with a KS significance higher than 0.33, and are also reasonably consistent 
with GRB970508, GRB971214, and GRB980703. Figure | shows where these three GRBs 
lie in comparison to the SFR model. 

We test whether number versus flux measurements from an instrument more sensitive 
than BATSE can differentiate the models. We consider an instrument with a flux cutoff of 
0.06 photons cm -2 s _1 , and perfect efficiency above this cutoff. Figure |2| shows that, while 
by construction the models predict identical numbers of bursts below the BATSE cutoff, 
they differ substantially at lower peak fluxes. Thus, the observed number of bursts per year 
can differentiate the models. As an interesting side note, in the SFR model, 48 percent of 
all GRBs in BATSE's field of view are seen with a flux cutoff of 0.42 photons cm -2 s -1 , and 
98 percent are above 0.06 photons cm -2 s~ x . 

We do a KS comparison of the CCD model to a variable number of GRBs chosen from 
the SFR distribution. Our results indicate that more sensitive instruments than BATSE 
may be able to differentiate our sample models using only the number-flux relation. One 
hundred and twenty-five bursts are enough for a two-sigma differentiation 95 percent of 
the time. The hypothetical instrument could make this observation in about 4 months, 
assuming 50 percent sky coverage and unit efficiency. 

While these results are encouraging, it seems unlikely that observations of flux alone 
would do this well in all cases. Based on N(< S), one model with many bright GRBs at 
high redshift might well be completely indistinguishable from another with dim GRBs at 
low redshift. In this case, redshift determinations would be necessary to differentiate the 
models. 

As shown in Figure ^|, the two models predict substantially different redshift 
distributions for the bursts. In order to determine the number of redshift measurements 
required to distinguish the two models, using BATSE's efficiency, we compared the 
CCD model to a variable number of bursts from the SFR distribution in a two-sided, 
two-dimensional KS test on the peak fluxes and redshifts. 

Our results indicate that about 25 redshifts are needed to distinguish the SFR and 
CCD models at the two sigma level 95 percent of the time. Our method is, however, a 
significant oversimplification: we assume that it is equally easy to measure any redshift. 
In practice, different measurement techniques result in an efficiency with strong and 
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Redshift 



Fig. 2. — The center panel shows points indicating the peak flux and redshift for 5000 
GRBs selected from the SFR distribution with the BATSE flux cutoff; the asterisks are 
GRB970508, GRB971214, and GRB980703. The left panel shows the number of GRBs per 
year in each peak flux range, divided into 20 evenly spaced logarithmic bins. The dot-dash 
vertical line shows the BATSE flux cutoff. The histogram in the bottom panel shows the 
expected number of bursts per year in each redshift range, divided into 30 evenly spaced 
logarithmic bins. The total number of bursts is generated by setting the number in each 
model equal to the number per year seen by BATSE, corrected for BATSE's effective coverage 
of 48.3 percent of the sky. In left and bottom panels, the solid histogram is the SFR model, 
and the dashed histogram is the CCD model. The dotted line in the left panel is the BATSE 
calatog. 
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complicated redshift dependence, with redshifts between ~ 1 and 2, where the bulk of star 
formation occurs, difficult to measure from the ground. 

Furthermore, two-dimensional KS tests are somewhat unreliable with small numbers of 
data points QPress et al. 1992j ). Parameters for the two-dimensional KS test are derived by 



numerical experiment, and assume that the number of samples is large enough for statistical 
treatment. Samples with fewer than about 20 data points give results that depend on the 
shape of the underlying distribution. Thus, our figures should be regarded only as lower 
bounds. Even so, we note that there are three measured redshifts out of fifteen accurate 
localizations. Given this efficiency, we would need more that 125 localizations to reach even 
this lower limit. 



5. Conclusion 

Our results indicate that the current data available from BATSE are insufficient to 
determine the redshift distribution of GRBs. The complicated mixing between luminosity 
and comoving redshift distributions make it impossible to significantly constrain either 
one at present. Our example models demonstrate that even a relatively small change in 
luminosity distribution can render both a star formation model and a constant comoving 
density model consistent with observations. 

The increased flux sensitivity provided by the next generation of satellites should be 
able to greatly reduce the parameter space available. If our comparison of the SFR and 
the CCD model is representative, more sensitive satellites may be able to confirm whether 
GRBs trace the star formation history even without localizing any bursts. However, redshift 
measurements provide considerably stronger constraints on GRB models. 

These results indicate that the best way to constrain the spatial distribution of GRBs 
is by directly measuring a sizable number of redshifts. As has been demonstrated by 
the Beppo-SAX localizations, this is best accomplished by rapid ground-based optical 
spectroscopic followup of accurately-positioned bursts, or optical spectroscopy of the host 
galaxies. A lower limit on the number of GRB localizations required to distinguish between 
disparate models, based on current efficiency, is more than a hundred. This is in excess of 
what will likely be achieved by continued operation of Beppo-SAX, and by the upcoming 
HETE-2 mission (which expects to localize <~ 50 bursts). It is likely, however, that more 
rapid notification to the ground of positions at the arc-minute or better level (~ 8' SAX 
positions currently take 8 hours before distribution, and refinement to 1' typically requires 
a day) will increase the efficiency for redshift determination by absorption spectroscopy of 
the optical afterglow. If this is the case, it is possible that the next generation of satellites 
may be able to test whether the spatial distribution of GRB sources follows star formation. 
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